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Earth-skimming UHE tau neutrinos have a chance to be detected by the Fluorescence Detector (FD) of Pierre 
Auger Observatory if their astrophysical flux is large enough. A detailed evaluation of the expected number of 
events is here performed for a wide class of neutrino flux models. 
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Neutrinos with energy above 10^^ eV are ex- 
pected to originate from the interaction of UHE 
cosmic rays with the Cosmic Microwave Back- 
ground (CMB) via the 7r-photoproduction, p + 
ICMB — * n-t-TT"*", the so-called cosmogenic neutri- 
nos The prediction for such a flux have been 
exhaustively discussed in several papers (see for 
example Ref.s |2I8| ^. 

High energy neutrinos are hardly detected, as 
they are almost completely shadowed by Earth 
and rarely interact with the atmosphere. In this 
framework, an interesting strategy for Vr detec- 
tion has been proposed in literature (see ^ for a 
complete list of references). For energy between 
10^* and 10^^ eV the r decay length is not much 
larger than the corresponding interaction range. 
Thus, an energetic r, produced by Charged Cur- 
rent (CC) i^T interaction not too deep under the 
surface of the Earth, has a chance to emerge in 
the atmosphere as an upgoing particle. Unlike 
r's, muons crossing the rock rapidly loose energy 
and decay. Almost horizontal i^r, just skimming 
the Earth surface, will cross an amount of rock 
of the order of their interaction length and thus 
will be able to produce a corresponding r, which 
might shower in the atmosphere and be detected. 

A detailed estimate of the number of possible 
upgoing T showers which the Fluorescence Detec- 
tor (FD) of Pierre Auger Observatory could de- 
tect is presented in Ref . ^ , where the predictions 
are analyzed with respect to their dependence on 



different neutrino fluxes and by using a new es- 
timate of neutrino-nucleon cross sections. In the 
present paper a brief review of the main results 
obtained in Ref. ^ is reported. 

Ultra High Energy protons, with energy above 
~ 10^° eV, travelling through the universe mostly 
loose their energy via the interaction with CMB 
radiation. The large amount of charged and 
neutral pions produced will eventually decay in 
charged leptons, neutrinos {cosmogenic neutrinos 
^), and high energy gamma rays. 




Figure 1. Cosmogenic neutrino fluxes as a func- 
tion of energy. 
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At the GeV energy range the extragalactic dif- 
fuse gamma-ray background was measured by the 
EGRET experiment. This measurement provides 
an upper bound for possible neutrino fluxes from 
pion production. In particular, it gives the ex- 
pected maximum flux of cosmogenic neutrinos 
from an initial spectrum of measured UHE pro- 
tons. It is worth noticing that, since at least part 
of UHECR are protons, the existence of cosmo- 
genic neutrinos is guaranteed, even if their flux is 
very uncertain. In Figure^the GZK neutrino flux 
for three possible scenarios is plotted. The thick 
solid line gives the case of an initial proton flux 
oc by assuming in addition that the EGRET 
flux is entirely due to 7r-photoproduction (GZK- 
H). The thin solid line shows the neutrino flux 
when the associated photons contribute only up 
to 20% in the EGRET flux (GZK-L). The dashed 
line stands for the conservative scenario of an ini- 
tial proton flux cx l/E"^ (GZK-WB). In this case 
the neutrino flux is compatible with the so-called 
Waxman-Bahcall limit. 
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Figure 2. Neutrino fluxes in exotic UHECR mod- 
els. 



Most of the models trying to explain highest 
energy cosmic rays {E > 10^° eV) in terms of 
exotic particles, predict a large associated flux of 
neutrinos. In Figurc|21 the expected neutrino flux 



for two of such scenarios is plotted. One of them 
is the model of new hadrons (NH), with mass 
M - 2 - 5 GeV, capable of generating UHECR 
events above GZK cutoff. In SUSY theories, for 
example, the new hadrons are bound states of 
light bottom squarks or gluinos and, once pro- 
duced in suitable astrophysical environments, can 
reach the Earth without signiflcant energy losses. 
In spite the production of new hadrons is a sub- 
dominant process, it generates a large number of 
neutrinos (see Figure 12). 

The dashed line in Figure |21 shows the neutrino 
flux for a Topological Defects model (TD) (for 
a review see IH]). In this case UHECR events 
with energy E > 10^*^ eV are explained in terms 
of 7's which are produced in the decay of heavy 
particles with mass of the order of io22~23 
in the previous case, the associated neutrino flux 
for this kind of models is extremely large. 

At energy above 1 GeV the interaction between 
neutrinos and the atoms of the rock is domi- 
nated by the process of Deep Inelastic Scatter- 
ing (DIS) on nucleons. Detectable leptons are 
produced through Charged Current interaction 
+ N + X. The total cross sec- 

tions can be written in terms of differential ones 
as follows 

o'-JhiEu) = " ^{E^,y)dy , (1) 

where E^, is the energy of the incoming neutrino, 
mi is the mass of the outgoing charged lepton and 
y is the inelasticity parameter, defined as 

ycc^i-§- , (2) 

with El the energy of the outgoing charged 
(for CC) or neutral (for NC) lepton. In the 
present analysis we have used the CTEQ6 [HI 
parton distribution functions in the DIS factor- 
ization scheme. The Q^-evolution is realized 
by the next-to-leading order Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi equations. 

Figure |2| shows a comparison between a 
CTEQ4-based parametrization of CC cross sec- 
tion and the corresponding calculation performed 
with CTEQ6. A substantial agreement is found 
up to 10^ GeV, whereas a discrepancy of at most 
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30% at IQi^ GeV is observed (CTEQ4 prediction 
being larger). In this range of energy the uncer- 
tainty due to the lack of knowledge of parton dis- 
tribution functions is expected to be overwhelm- 
ing. 
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Figure 3. The total i^T-nucleon CC cross sections 
based on both CTEQ4 |Z| and CTEQ6 are 
reported. 



Following the formalism developed in Ref. [H], 
let be an isotropic flux of + Vr- The differ- 
ential flux of charged leptons emerging from the 
Earth surface with energy E-r is given by 



d<^r{Er,^) 

dErdVL 



dE, 



d^,{E,,n) 
dE^dn 
K{E„ n; Er) , 



(3) 



where K{E^, 0; E^) is the probability that an in- 
coming neutrino crossing the Earth with energy 
En and nadir angle 9 produces a lepton emerging 
with energy E^ (see Figure QJ. In Eq.©, due to 
the very high energy of v^, we can assume that 
in the process Vr + N ^ t + X the charged 
lepton is produced along the neutrino direction. 




Figure 4. A neutrino Vr crosses the Earth with 
energy E^ at a nadir angle 9 and azimuth angle 
(j). Then it travels a distance z before converting 
into a charged lepton r, which emerges from the 
Earth surface with energy Ej- 



This process can occur if and only if the following 
conditions are fulfilled: 

a) the Vr with energy E^ has to survive along a 
distance z through the Earth; 

b) the neutrino converts into a r in the interval 

z,z + dz; 

c) the created lepton emerges from the Earth 
before decaying. 

By using the expressions for all the above prob- 
abilities one obtains 



K{E,, 9; Er) = 



X exp 



Er{f3r+lrEr] 
1 1 



{F(E,,Er)f 



CTrPrQs \Er E^{E„ 

2Rq COS 9 a''cc{Eu)QsN a] , 



(4) 



where ^ = uj + a'^'^iE^) Na/ Pr- As exten- 
sively discussed in Ref. @] the parameters (3r — 
0.71-10-6 cm^ g-i and 7^ ~ 0.35-10-1* cm^ g-i 
GeV-i fairly describe the r energy loss in matter. 
Moreover, the tau lepton, once produced, carries 
an average energy E^{En) = {1- < ycc >)Ey. 

Eq. (0J) leads to the total rate of upgoing t's 
showering on the Auger detector, and thus poten- 
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tially detectable by the FD: 



= 2nSD / dE^ / dEr 

dt Je""" J e*^ 



X 1 — cxp 



dE^dn 
H nir 



K(E,, 0; Er) 



CTr Et 



e co8 6d{cos9) , (5) 



where we have used the isotropy of the consid- 
ered neutrino flux. In Eq.(|3J) the quantity 5 — 
3000 km^ is the geometrical area covered by the 
Auger apparatus, D ^ 10% is the duty cycle for 
fluorescence detection, ~ 10"'^* eV is the en- 
ergy threshold for the fluorescence process, and 
jTjmin jg ^Yie minimum neutrino energy capable of 
producing a r at detection threshold. The quan- 
tity E^"'^ is the endpoint of the neutrino flux. 

The exponential term in the r.h.s. of Eq.(|Sl) ac- 
counts for the decay probability of a t (showering 
probability) in a distance H from the emerging 
point on the Earth surface. 

In Eq. (|SJ) the integration over cos 9 can be eas- 
ily performed and this yields to the definition of 
the effective aperture of the apparatus, A{E^), as 



dNr 
"df 



= D 



In Figure [S] the quantity A{Ei,) is plotted versus 
the neutrino energy for different values of E^^ and 
H. 

The number of r-shower Earth-skimming 
events expected per year at the FD detector for 
the different neutrino fluxes result to be: 0.02 
(GZK-WB), 0.04 (GZK-L), 0.09 (GZK-H), 0.11 
(TD), 0.25 (NH). 

As well known, neutrino induced extensive air 
showers (see Ref. |3 for an extended discussion) 
can be disentangled by the ordinary cosmic ray 
background only for very inclined and deep show- 
ers. An estimate of the expected number of down- 
going events within 30 km from the FD detector 
results to be comparable with the upgoing one 
only for zenith angles larger than 70°. 
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Figure 5. The effective aperture A{E,y) is here 
plotted for E^J" = lO^"^ eV and H ^ 30 km, and 
for El'' = 10^8 eV with H = 20, 30, and 40 km, 
respectively. 
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